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SUMMARY 

This paper describes a new geometric analysis procedure for wing sections. This procedure is based 
on the normal mode analysis for continuous functions. A set of special shape functions is introduced to 
represent the geometry of the wing section. The generators of the NACA 4-digit airfoils were included 
in this set of shape functions. It is found that the supercritical wing section, Korn airfoil, could be well 
represented by a set of ten shape functions. Preliminary results showed that the number of parameters 
to define a wing section could be greatly reduced to about ten. Hence, the present research clearly 
advances the airfoil design technology by reducing the number of design variables. 

STATEMENT OF PROBLEM 

In the optimization procedure, wing sections are perturbed by linear combinations of shape functions 
to render an optimal geometry which satisfies certain mission requirements. The final optimal geometry 
and the efficiency of the design tool depends upon the choice and the number of shape functions used. 
The selection of shape functions to suitably represent wing sections is one of the fundamental problems 
to airfoil designers. 

Mathematically, any continuous function defined on a closed interval can be represented by an 
infinite series of normal modes which form a complete set of bases. The set of Fourier sine functions 
is an example of such a complete set. However, practical application does not call for an infinite series 
but a finite partial sum to approximate the geometry of the wing section within a prescribed tolerance 
in some norm. 

There are several well known shape functions for wing section modification introduced in recent 
papers. Hicks functions, Wagner functions, and aerofunctions are some of them (refs. 1-3). Although 
these shape functions are useful, there are no systematic guidelines to effectively select these shape 
functions. The number of shape functions used in current design optimizations is easily more than 
20 — too many to be efficient. 

The present research pondered this problem and proposed a solution to it via normal mode analysis 
(ref. 4). 
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NORMAL MODE ANALYSIS 


In the following, we briefly describe the normal mode analysis of the continuous function. 

An inner product of two integrable functions, g and h, on the closed unit interval is denoted as 
(g,h) and is defined by the relation 


l 

(g,h) = J y{x) h(x)dx 
0 

A norm of / is defined as the square root of (/,/). Two functions, g and h, are orthogonal to each 
other if ( g,h ) = 0. A set of functions {fi{x)} is orthogonal if (f z , fj) = 0 for i / j and ^ 0 

for all i. A set of functions {f z (x)} is orthonormal if it is orthogonal and the length of each member 
is unity, ( fi,f z ) = 1. A set of functions is linearly independent if the norm of any linear combination 
of its members is zero, provided that the coefficient of each term of the combination is zero. 

Let {< 7 j (x)} be a set of functions which is linearly independent. The Gram-Schmidt orthonormal- 
ization process for the set {<7j(x)} can be described as follows. 

First, an orthogonal set {/i;(x)} is formed from the following relations: 

hi(x) = g\(x) 
h 2 {x) = g 2 {x) - a 2 \hi{x) 


71 — 1 

h n (x) = g n (x) - a 7 ii f H(x) 


7 — 1 


where a ni = (g n , h x )/(/;», hj) is the projection of g n in the direction of /»,. Finally, the orthonormal set 
{fi(x)} is found by normalizing /<j(x) as follows: 

fi{x) = hi(x)/(hi,hi) 


The least-squares problem is to seek a vector f(x) in the subspace spanned by {.g,(x)J which lies 
the shortest distance from the given vector f(x). Namely, the norm of error function, / - /, is as small 
as possible. This least-squares problem has a unique solution: 

/(*) = Y. b ifi( x ) 

7 — 1 

where the orthonormal coefficient is 6, = (/,/,). The function / is called the reconstruction of function 
/ via the normal mode analysis. 





PRACTICAL APPLICATION 


An airfoil geometry is usually defined by a discrete set of points. There are three steps involved 
in the process to normalize the airfoil geometry into the canonical form which requires that the leading 
edge be located at the origin of the coordinate system and the trailing edge at unit length of the horizontal 
axis. First, for an airfoil with finite thickness at the trailing edge, some work has to be done to close the 
trailing edge of the airfoil. Second, the leading edge of the airfoil has to be determined. The leading 
edge of the airfoil is determined as the farthest point from the trailing edge which lies on the circle 
passing through the three leading points of the discrete set which are farthest away from the trailing 
edge. Finally, a scaling and a rigid body transformation (rotation and translation) may be needed to 
map the airfoil geometry onto the canonical form. 

A complete geometry of wing section entails two parts, upper and lower surfaces. They can be 
combined to define the camber and thickness distributions by the following relations. Namely, let the 
net {x^} be a discrete set of abscissas on the closed unit interval and {uj} and {t ! j} the upper and lower 
surface functions defined on the net. Then, the camber distribution {cj} is given by the relation 

Cj = 0.5(itj + v t ) 

and the thickness distribution {/;} by the relation 

U = Uj - Vi 

Conventional Airfoil 

There are several books which include airfoil definitions in the literature (refs. 5-7). The NACA 
0012 airfoil is a typical example for conventional airfoils. It falls in the category of NACA 4-digit 
airfoils (ref. 5) which are defined by the following equation: 

y = -^(0.2969Vx - 0.1 26x - 0.3516x 2 + 0.2843x 3 - 0.1015x 4 ) 

where t, is the maximum thickness expressed as a fraction of the chord which is of unit length. It is 
clear that function y can be rewritten in terms of the following four shape functions: 


II 

1 

H 

(1) 

g 2 {x) = x(l - x) 

(2) 

93 (x) = x 2 (l - x) 

(3) 

34 (x) = z 3 (l - x) 

(4) 


where x is on the closed unit interval [0,1]. These four functions are referred to as the generators of 
the NACA 4-digit airfoils. 

The geometry of the original NACA 0012 airfoil defined above is not closed at the trailing edge. 
An effort was made to close the trailing edge of the airfoil. Figure 1 shows the contour of the normalized 
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NACA 0012 airfoil. The set of shape functions, {</j(x),z = 1,4}, is plotted in figure 2(a). The set 
of orthonormalized duals, {f t (x),i — 1,4}, of the shape functions is displayed in figure 2(b). The 
orthonormal coefficients of the upper surface of the NACA 0012 with respect to the orthonormal duals, 
{fi(x),i = 1,4}, are shown in figure 3. The guppy plot of the reconstructed airfoil, generated by 
the orthonormal coefficients and the orthonormal duals, is displayed along with the input airfoil for 
comparison in figure 4. The maximum norm of the error function of the NACA 0012 airfoil was about 
0. 18E-4. This showed that the NACA 0012 airfoil was practically recovered through the normal mode 
analysis. 


Supercritical Airfoil 

The extension of the NACA shape functions to represent a supercritical wing section is now 
presented. The following six shape functions which are in the homotopic neighborhood of the NACA 
shape functions were added to the NACA listing. They are as follows: 


g 5 {x) = x 4 (l - x) 

(5) 

g e {x) - x 5 (l - x) 

(6) 

g 7 {x) = tyx - y/x 

(7) 

g 8 (x) = $[x - yfx 

(8) 

g 9 (x) = ific - tyx 

(9) 

9lo(x) = yfx - y/x 

(10) 


Now there are ten shape functions in this new NACA collection. 

The supercritical wing section, Korn airfoil, was chosen to be studied by these ten shape functions. 
The Korn airfoil was not defined by an algebraic equation but was devised by an inverse design method 
on the hodograph plane. The geometry definition of the Korn airfoil can be found in reference 6. Its 
trailing edge is closed. The leading edge of the airfoil was not given in the original discrete set for the 
airfoil geometry. The leading edge is determined to be (—0.000053,0.002228) and the pitching angle 
between the leading edge and trailing edge was calculated as 0.127663 degree (nose up). Figure 5 
shows the geometry of the normalized Korn airfoil. The set of shape functions, {gi(x),i = 5,10}, 
is plotted in figure 6(a). The set of orthonormal modes, {fj(x),i = 5,10}, of the shape functions is 
shown in figure 6(b). The orthonormal coefficients of the upper surface and lower surface of the Korn 
airfoil corresponding to the orthonormal modes, {fi{x),i = 1, 10}, are shown in figure 7. The guppy 
plots of the reconstructed Korn airfoil along with the input Kom airfoil are displayed in figure 8. The 
maximum norm of the error function of the Korn airfoil was about 5.70E-4. It shows that the Korn 
airfoil can practically be represented by these ten shape functions. 

CONCLUDING REMARKS 


A geometric analysis was described for practical wing sections. The approach of the analysis was 
based on the approximation theory. The least-squares approximant to the geometry of the wing section 
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was determined via the normal mode analysis. A suitable set of shape functions was newly gathered 
to generate a set of orthonormal modes through the Gram-Schmidt orthonormalization process. The 
generators of the NACA 4-digit airfoils were included in this set of NACA bumps. The geometry 
contour of the NACA 0012 airfoil was represented by a combination of only four normal modes within 
a very tight tolerance, whereas the Korn airfoil, a supercritical wing section, was well represented by the 
set of only ten orthonormal modes. Preliminary results showed that the number of parameters to define 
a practical wing section can be reduced to about ten of these newly identified shape functions. The 
significance of the present research is to provide a promising method for advancing the wing section 
design in reducing the number of design variables. 

Additional work remains to be done. One area of particular interest is to apply the present geometric 
analysis to other advanced airfoils to further confirm the usefulness of the NACA bumps. Another is to 
make direct comparison of airfoils with their reconstructed duals in terms of aerodynamic characteristics. 
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irfoil with the trailing edge closed. 







Coefficient 



Figure 3. The normal mode coefficients of the upper surface of the NACA 0012 airfoil corresponding 
to the normal modes, {/ 7 (.r), i = 1,4}. 
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Input/Reconstructed NACA 0012 Airfoils 


Input Airfoil 
Reconstructed Airfoil 



Figure 4. The guppy plots of the input NACA 0012 airfoil and its reconstructed airfoil via the normal 
mode analysis. 
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il with the fitted-in leading edge. 








NACA Normal Mode Analysis for Kom Airfoil 
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Figure 7. The normal mode coefficients of the upper and the lower surfaces of the Korn airfoil corre- 
sponding to the normal modes, {/,(./'), i = 1,10}. 
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